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electron- and hole-like quasi-particles) on the superconducting layers and investigate its influence 
on the quasi-particle current between the layers. This nonequilibrium effect leads to shifts in the 
current-voltage curves of the stack. The theory is applied to the interpretation of recent current 
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I. INTRODUCTION 

In the strongly anisotropic cuprate superconductors 
like Bi2Sr2CaCu20s+5 (BSCCO) the superconducting 
CuC>2 layers together with the intermediate insulating 
material form a stack of Josephson junctions. In the 
presence of a bias current perpendicular to the layers 
the intrinsic Josephson effect manifests itself in a multi- 
branch structure of the IV-curves^2i2iii^ where on each 
branch a different number of junctions is in the resistive 
state. 

In the resistive state high frequency Josephson oscil- 
lations at a finite dc-voltage are accompanied by a dc- 
current which is carried mainly by unpaired electrons 
(quasi-particles), while in the superconducting state it 
is carried mainly by Cooper pairs. Therefore, on a su- 
perconducting layer between a resistive and a supercon- 
ducting junction the bias current has to change its char- 
acter from quasi-particle current to supercurrent. This 
creates a non-equilibrium state on this layer with a fi- 
nite quasi-particle charge and a change of the condensate 
charge. The quasi-particle charge is characterised by a 
distribution function with different numbers of electron- 
like and holelike quasi-particles, while the change of the 
condensate charge, the number of paired electrons, can 
be described by a shift of the chemical potential of the 
condensate. 

In recent experiments such non-equilibrium effects 
have been observed in layered d-wave superconductors^ 
In a first type of experiments Shapiro steps produced by 
high-frequency irradiation have been measured in mesa 
structures of BSCCO with gold contacts. Here a shift 
of the step voltage from its canonical value hf/ (2e) was 
observed, which can be traced back to a change of the 
contact resistance due to quasi-particle charge on the 
first superconducting layer. In another type of experi- 
ments current-voltage curves have been investigated for 
two mesas structured close to each other on the same 
base crystal (see Fig. 7 below). Here an influence of the 
current through one mesa on the voltage drop through 
the other mesa has been measured which is caused by 
charge imbalance on the first common superconducting 
layer of the base crystal. In a recent paper— we have ex- 



plained these effects by using a semi-phenomenological 
approach based on a microscopic non-equilibrium Green 
function theory 8 ^ for layered superconductors. 

In this paper we will present the full microscopic the- 
ory for stationary charge-imbalance effects in intrinsic 
Josephson systems. In order to be specific, we apply 
the theory first to an experiment, which in this form has 
not yet been done, but is conceptual simpler then the 
double-mesa injection experiment mentioned above. We 
consider a mesa as shown in Fig. 1 with two normal elec- 
trodes on top. Through one electrode a stationary bias 
current is applied creating a charge imbalance on the first 
superconducting layer. At the other electrode a voltage 
is measured as function of the bias current through the 
first electrode with zero current through the second elec- 
trode. This experiment is similar to the classical exper- 
iment by Clarke ^fi where a strong current injected into 
a bulk superconductor creates quasi-particles. By two 
other electrodes, one normal junction and one Joseph- 
son junction, a voltage is detected at zero current which 
measures the difference between the chemical potential of 
quasi-particles and condensate. In our case the junction 
between the normal electrode and the first superconduct- 
ing layer is the normal junction, the coupling between the 
first superconducting layer and the next superconducting 
layers is the Josephson junction. This experiment allows 
to determine the charge imbalance and (with further the- 
oretical input) its relaxation rate. After this investigation 
we return to the discussion of the double-mesa injection 
experiment. 

Nonequilibrium effects in bulk- superconductors have 
been studied a long time ago. For an overview see 
the text-book by Tinkhamii and the review-articles 
edited by Langenberg and Larking The basic concepts 
of charge-imbalance have been developed by Tinkham 
and Clarke; 10 i 13 Pethick and Smith p11 The theory has 
been worked out using a non-equilibrium Green func- 
tion approach by Schmid and Sch6n^£*i& and Larkin 
and Ovchinnikoviii^ Microscopic nonequilibrium the- 
ory of layered superconductors was considered first by 
Artemenkoi^ and later by Graf et ali^S 

Nonequilibrium effects in intrinsic Josephson sys- 
tem have been investigated theoretically already 
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FIG. 1: Mesa-structure used for 4-point measurements con- 
sisting of a stack of superconducting layers with two normal 
electrodes on top. 



in various contexts using different met hods and 
approximations, 21 . 22 ! 23 ! 24 ! 25 . 26 ! 27 ! 28 ! 29 In Refs. I2H2I the 
influence of charge-coupling on Josephson plasma os- 
cillations has been investigated, in the approach used 
in RefsiSiSiSi a systematic perturbation theory in the 
scalar potential on the layers has been performed. In 
these theories charge imbalance is considered only indi- 
rectly as far as it is induced by fluctuations of the scalar 
potential. The present theory which uses non-equilibrium 
Green functions is more general. Here charge imbalance 
is taken into account as an independent degree of free- 
dom, and therefore the results are different from those of 
earlier treatments. 

Here we apply this theory to stationary processes in 
layered d-wave superconductors. Due to the weak cou- 
pling between the superconducting layers the total charge 
created by the bias-current is non-zero. Therefore we 
have to treat charge fluctuations of the condensate and 
quasi-particles as independent quantities. The d-wave 
character of the order-parameter allows quasi-particlc 
tunnelling also at low voltages and temperatures. Fur- 
thermore it allows the relaxation of charge-imbalance by 
elastic scattering processes (impurity scattering within a 
layer). On the other hand the weak coupling between the 
layers leads to a partial decoupling of the kinetic equa- 
tions for the Green functions in different layers, making 
this system a particularly simple example for the appli- 
cation of nonequilibrium theory in superconductors. 

In the following section we start with general consider- 
ations concerning the Josephson effect in layered super- 
conductors and the definition of charge-imbalance. We 
then formulate the basic kinetic equations describing sta- 
tionary nonequilibrium effects in intrinsic Josephson sys- 
tems and calculate the charge imbalance microscopically. 
We show how these effects can be measured by discussing 
the above mentioned 4-point experiment and the cou- 
pled two-mesa system. Finally we discuss the influence 
of charge-imbalance on the IV-curves of a stack of junc- 
tions. 



II. GENERALISED JOSEPHSON EQUATIONS 

We consider a system of superconducting layers with 
indices n and order parameter A„(t) = | A| exp(ix„(i)) 
with time-dependent phase Xn{t)- We define the gauge 
invariant phase difference as: 



7n,n+l(*) = Xn(t) - Xn+l(t) - 



2c 



n+1 



dzA z (z,t), (1) 



where A z (z, t) is the vector potential in the barrier. Here 
e denotes the elementary charge. The charge of the elec- 
tron is — e. 

For the time derivative of 7„ in +i we obtain the gener- 
alized Josephson relation: 



^7n,raH 

dt 



^(v 
h V n 



„+l + $„+i - 



Here 



V, 



n,n+l 



n+1 



dzE z (z, t) 



(2) 



(3) 



is the voltage and <fr„(i) is the so-called gauge invariant 
scalar potential defined by 



<S>n(t)=<Pn(t)-—Xn{t), 

Ze 



(4) 



where <j> n (t) is the electrical scalar potential. 

The quantity /ry„, n+ i = 2e(V n , n +i + $n+l - $n) is the 
total energy required to transfer a Cooper pair from layer 
n to n+1, e$„ can be considered as the shift of the chem- 
ical potential of the superconducting condensate with re- 
spect to an average chemical potential /i, i.e. the number 
of particles in the condensate is controlled by /i+e$ n . For 
equilibrium superconductors %Xn — 2e</>„, <!>„ = 0, and 
one has the usual Josephson relation ^7 n>n +i = 2eV n _ n +i. 

The generalised Josephson relation J2J) for a stack of 
Josephson junctions can be written in another way show- 
ing explicitly the electric field distribution between the 
layers. We split the total charge fluctuation on a layer 
5p n = + (5/?^ into a contribution bp c n from particles 
in the condensate and a contribution <5p* from unpaired 
electrons (quasi-particles) describing charge-imbalance. 
The charge fluctuation of the condensate can be ex- 
pressed directly by the change of the chemical potential 
on layer n: 



6p c n = -2e 2 7V(0)$„, 



(5) 



where N(0) is the (two-dimensional) density of states for 
one spin direction at the Fermi energy. It is convenient 
to express also the fluctuation of quasi-particle charge by 
some quasi-particle potential ty n , defining 



5p q n = 2e 2 iV(0)*„. 



(6) 



Then we obtain for the total charge density fluctuation: 
6p n = -2e 2 AT(0)($ n -* n ). (7) 
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While in bulk superconductors charge neutrality leads 
to $ = 'J, this is not the case in weakly coupled layered 
superconductors. 

With help of Q and the Maxwell equation (d is the 
distance between the layers) 

Sp n = —T-(Vn,n+l - V n -\,n) (8) 

the generalized Josephson relation now reads: 

h . 

7^7n,n+l = (1 + 2a)V n .n+l 

2e 

- a(V n -i tn + V n+ l, n+2 ) + *n+l - *n, (9) 

with a — eeo/ (2e 2 N(0)d). It shows that the Joseph- 
son oscillation frequency is determined not only by the 
voltage in the same junction but also by the voltages 
in neighboring junctions. The coefficient a has already 
been introduced by Koyama and Tachiki(2i however in 
their theory charge imbalance effects described by the 
quasi-particle potential ty n have not been considered. 

For a barrier in the superconducting state the time 
average (7,1,71+1) vanishes, while for a barrier in the re- 
sistive state ft(7„ iIl+ i)/2 is the electrochemical potential 
difference. A similar equation also holds for the contact 
with the normal electrode. Denoting the normal elec- 
trode, which is assumed to be in equilibrium by n = 0, 
we have 

^-7o,i = Vb,i + $i = (1 + a)V ,i + aVi,2 + (10) 

In the case of a stack in contact with a normal electrode 
on top, where all internal barriers are in the supercon- 
ducting state we find for the first superconducting layer 
$!~oy + *i. This result follows from Eqs. m 
the limit a<l. Note that the total voltage of the stack 
is V = E n 7n,n+lV(2e). 

In order to describe transport between the layers we 
need also an expression for the current density. In our 
previous papers^! we have used 

jn,n+l = jc sin 7„, n+ i 

+ ^ ±i (^7n,„ + l + *n-*n +1 ), (H) 

which is approximately valid in the stationary case (no 
displacement current). The first term is the current den- 
sity of Cooper pairs. The rest is the quasi-particle current 
density, which is driven not only by the electrochemical 
potential difference, but contains a diffusion term pro- 
portional to the quasi-particle density difference. In the 
following the current expression will be derived from the 
microscopic theory. 

Finally we need some theory describing the creation 
and relaxation of quasi-particle charge in order to calcu- 
late the charge-imbalance potential In our previous 
papers we have used a relaxation time approximation. In 
the stationary case we obtained 

*n - r q { 3 l_ 1;n - ^ >n+1 )/[2e 2 iV(Q)] (12) 



describing the balance between charge-imbalance cre- 
ation by the in- and outgoing quasi-particle currents and 
relaxation into thermal equilibrium inside the layer. It 
will be the main task of this paper to derive such a rela- 
tion from a microscopic theory. 



III. KINETIC EQUATIONS 

For the microscopic description we start from equa- 
tions of motion for the spectral ( retarded and ad- 
vanced) and Keldysh Green functions in Nambu space 
(for details see the appendix^ . Nonstationary cor- 
rections to the spectral functions are small, and these 
functions can be taken in equilibrium. The Keldysh 
function contains the necessary information about the 
nonequilibrium distribution functions. In our model we 
consider d-wave superconductivity on each layer within 
the BCS-approximation, elastic impurity scattering in- 
side the layer and tunnelling between neighboring layers 
in lowest order. Then it is possible to formulate equa- 
tions of motion for the Green function G^{ji,tx,ti) for 
each layer n with a tunnelling self-energy containing the 
coupling to the neighboring layers. 

Furthermore we introduce quasiclassicle Green func- 
tions, which are obtained by integrating the Green func- 
tion over the kinetic energy — ek—fi- Performing in ad- 
dition a Fourier transformation with respect to the time 
difference r = t\ — t-x at fixed central time t= [t\ + t2)/2 
we obtain (see the appendix[5J) the following kinetic equa- 
tion for the quasiclassical Keldysh function g„ (k, t, e) in 
Nambu space: 

(ll + e ) T3 ^ + (^- £ )^ T3 
- {h(k,t)g K - g K h{k,t)} 
= {a R g K + a K g A -g Ii a K -g K a A } (13) 

depending on the central time t, the direction k of the 
momentum, and the frequency e. Here 

h(k,t) = -e$ n (t) - A(k), 

**> " U> f ) ■ <"> 

where & n (t) is the gauge invariant scalar potential 
Eq.(@J. A(/c) is the superconducting order parameter 
with d-wave symmetry. We have applied a gauge trans- 
formation such that on each layer the superconducting 
order parameter has the same phase (in our case it is 
real). The quantities a are self-energies due to impurity 
scattering and tunnelling, which will be specified later. 
The curly brackets denote a convolution in time and fre- 
quency space. In the stationary case, which we consider 
in the following, we will neglect the time dependence. 
Then the convolutions are simple products, the depen- 
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dence on drops out, and we obtain the basic equa- 

tion 



= [a R g K + a K g A -g R a K -g K a A ]. (15) 



This is actually an equation for the nonequilibrium distri- 
bution functions contained in the Keldysh Green function 
g K and the Keldysh selfenergy <jk- 

In the following we will use the ansatz by Larkin and 
O vchinnikov 



with 



g K ={g R d-dg A ) 



d = [3 + af 3 



(16) 



containing two distribution functions /3 and a. In the 
stationary nonequilibrium case we can neglect the con- 
volution and replace it by products: 



fJ 



K 



(g R -g A )P + (g R T 3 -T 3 g A )a, 



(17) 



where for the retarded and advanced functions g R ' (e) 
the equilibrium functions can be used. In equilibrium 



/3(e) = tanh(e/(2T)), a(e) = 0. 



(18) 



Returning now to the equation of motion Ea. l|15|) 
and using the traditional notation (see appendix A) for 
matrices in Nambu space, gu —: g,g\2 ='■ /, <?2i =: 
— / + ,<?22 =: g, the stationary part of the equations of 
motion for the normal and anomalous Keldysh Green 
functions read explicitly 





-A(f +K 


-f K ) 


= hi, 


(19) 
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-A(/ +A ' 


-f K ) 


= I22, 


(20) 


2ef K 


+ Mg K 


-g K ) 


= h2, 


(21) 


2ef +K 


+ Mg K 


-g K ) 


= hl, 


(22) 



with the abbreviation for the tunnelling and scattering 
integrals: 



T r A-R ?,K , ~K -A ~R~K 

lap = [<y g g -go- 



g K e A ~ 



a/3' 



As we shell see later (/) = (au) = (021) = for d-wave 
symmetry, and we obtain: 



hi 


= <Tll9 


"5^4 




9 K °ti 


I22 


^ R 7, K 

= G 22g 


-g R °^ 


K -A 

- <?22g - 


-K A 
g °22 


I12 


= Put 


-J a 22- 


L ^ K f A 

^11/ - 


- J a 22 


hi 


= -°22i 


f K _j_ f +R 


~K K 
a l\ ~ °22 


f +A + . 



The scattering terms on the r.h.s have different mean- 
ing: those which are proportional to the function itself 
can be considered as self-energy. Therefore we add these 
terms to the l.h.s. of Eas. ltzTT|) and J22J), writing: 



2ef K + A(g K - g K ) = J v 



(24) 



with 



Rea^ 



R,A 



(25) 



R,A 



where we used the equilibrium values for <7 n ' 
and u A x = {p\x)*- Then the correction to e is just the 
real part of the scattering self-energy evaluated below. 
The scattering term on the r.h.s. is then 



Jv 



fRK , K t A 



T f+ R ^ K ^ f 

•'21 = J cr u - <7 2 2J 



K r + A 



(26) 



Now we substitute the anomalous Keldysh functions 
f K and f +K in Eqs.JHI HUJ) and obtain 



hi 



•hi) = 



hi 



'22 







(27) 



(28) 



These equations will be used in the following to deter- 
mine the distribution functions a(e) and /3(e) describ- 
ing the non-equilibrium state. In particular, Eg. 1271) 
describes the balance between the relaxation of quasi- 
particle charge due to impurity scattering within a layer 
and its creation by the tunnelling current (see Ea. 143(1 
for the result in a special case). 



IV. CALCULATION OF THE CHARGE 
IMBALANCE 

We consider potential scattering from randomly dis- 
tributed impurity centers and tunnelling to the neigh- 
boring layers in lowest order. In order to calculate the 
self-energies and scattering integrals we need expressions 
for the retarded and advanced Green functions. Here we 
use the results from equilibrium (this can be justified by 
studying the corresponding equations of motion in the 
low frequency limit). 



A. Retarded and advanced Green functions 

The retarded (advanced) quasi-classical Green func- 
tions in equilibrium have the form: 

g R (ke)= , e + llt g A = -{g R )\ 



f R (k,e) 



{e + i lt f - A 2 (k) 

AW 

{e + i lt f - A 2 (fc)' 



r = -(/T, (29) 



with — i~f e — afi{e). Here we use the sign-convention 
Imyf > 0. Furthermore we have for the equilibrium 
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functions f+ R > A = f R > A , g R ' A = -g R > A . In the following 
we will also need the combinations: 

u(k,e) = ~(g R (k,e)-g A (k,e)), 

v(k,e) = ±(f R (k,e)-f A (k,e)) 7 

w(k,e) = ~(f R (k,e) + f A (k,e)). (30) 

The (even) spectral function u(k, e) is the tunnelling den- 
sity of states. With this notation we obtain for the 
Keldysh functions Eq. i|17|) 

g K (k 1 e) = 2u(k,e)(a(k 1 e)+f3(k 1 e)), 

g K (k, e) = 2u(k, e)(-f3(k, e) + a(k, e)). (31) 

B. Potential scattering 

The self-energies in Born approximation are given by: 

= Wo\ 2 G R ^ K (k,e) = -i Vp{ g R ^ { e)) 

k' 

(32) 

with v p = c7rA^(0)|Vo| 2 . Here Vq is the scattering poten- 
tial and c the impurity concentration. For comparison in 
the numerical results we also use the t-matrix approxi- 
mation in the strong scattering limit. This is obtained 
by replacing v p by v p /\g R (e)\ 2 . 

For the calculation of the Keldysh component of the 
self-energy we use the ansatz by Larkin and Ovchinnikov 
Eq. (|T7|) and obtain 

a* (e) = -2iu p (u{e)l3(e))T 3 - 2iu p (v{e)a(e))T Q . (33) 

Using these results we obtain for the potential part of 
the scattering terms: 

III = -I P 22 = Aw p [u{e){u{e)a{e)) - (u(e))u(e)a(e)] , 
J i2 ~ J 2i = 8ivpv(e){u(e)a(e,t)). (34) 

C. Tunneling 

We describe the coupling between neighboring layers 
by the usual tunnelling Hamiltonian 

h= E (^ fc ^* /:+1 ^ (2 ' t) ci +1 , fc v^+^.c.), 

ncrkk' 

(35) 

which after the gauge transformation becomes 

H = (W^»' (t) 4 +1)fe , CTCn , fcCT + h.c). (36) 

ncrkk' 

Then we obtain for the contribution to the tunnelling self- 
energy from the coupling between layers n and n' = nil: 



*t(e, t) n ,n> = / d(h - t 2 y^~^ \tk'k\ 2 

^ k> (37) 

e- ir37 "'"' (tl)/2 G(fc',ti,t 2 )e +lT37 "^' ( * 2)/2 . 

If we neglect additional phase fluctuation, the phase dif- 
ference is J n ,n'(t) = 7°,n' + ^n,n'*> where 7° „, is the 
static phase difference, which is determined by the super- 
current, and fi n , n ' is proportional to the electrochemical 
potential difference between the two layers. 

For the tunnelling matrix element we have to as- 
sume a partial conservation of momentum in order to 
obtain Josephson coupling between different supercon- 
ducting layers with d-wave order parameter. This, 
however, is not relevant for the tunnelling of quasi- 
particles, which will be discussed here primarily. In 
the following we will model the tunnelling matrixelement 
near the Fermi surface by irN(0)\tk>k\ 2 = v t (k : k ) with 
{A(k)u t (k,k')A(k')} ^0. 

The most simple case to study nonequilibrium effects 
is the tunnelling between a normal electrode and a super- 
conducting layer, which will be discussed in the following. 
Results for quasi-particle tunnelling between two super- 
conducting layers will be given in the appendix. In the 
case of tunnelling with a normal elctrode only the normal 
Green functions contribute to the tunnelling self-energy, 
which is time- independent. Writing v = Sl„.„//2 we ob- 
tain 

<xf(M(n,«') =-\(iyt9 R {e-v))' -{v t g R {e + v))')T 
-^t 9 R (e ~ v))' + (u t g R (e + v))')r s , (38) 

<rf(fc,e) („,„,, = (39) 

- i{[v t u'{e - v)(3'(e - v))' - {v t u'(e + v)(3'(e + v))'^t 

- i([v t u'(e - v)p'(e - v))' + {u t u'(e + v)(3'(e + v))')t 3 . 

The prime denotes spectral functions and distribution 
functions on layer n' and {vtg R )' = (vt(k, k')g R (k/ , e))k> 
denotes an average over the direction of k' . 

These results will be used to calculate the tunnelling 
contributions to the scattering terms I|27I28[1 in the ki- 
netic equation. Denoting the superconducting layer by 
n, the normal layer by n' , replacing the spectral function 
on the normal layer by u'(e) = 1, and neglecting terms 
of order vta we find 

/*? - = 4i W(e -v)- P'{e + v))u(e)Y, 

III + I22 1 = 4*M/3'(e -v) + p'(e + v) - (3{e)u{e))' , 

J% - JTi = 4i{v t v(e)(p'(e - v) - (3'(e + «))>'. (40) 

In order to determine the charge-imbalance on the su- 
perconducting layer in the stack we also need the quasi- 
particle contribution from tunnelling into the neighboring 
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superconducting layer with the barrier in the supercon- 
ducting state. As will be shown in the Appendix, this 
contribution vanishes if we neglect terms of order v t a. 



D. Solution of the kinetic equation 

Inserting now the different contributions to the scatter- 
ing term from potential scattering Eq. (|34|l and tunnelling 
Eq. (|40|l into Eq. i|27|) we obtain the following equation de- 
termining the distribution function a(k, e) on the super- 
conducting layer: 



2p p [u(e)(u(e)a(e)) - (u(e))u(e)a(e)] 
Av(e) 



(41) 



2v„ 



u(e)a(e)) 



= -(u(e) - jv(e))(u t (k, k') [(3>(e - v) - f3'(e + »))] )'. 

The r.h.s. of this equation describes charge imbalance 
generation by tunnelling of quasi particles. The l.h.s. 
describes charge imbalance relaxation due to impurity 
scattering. Taking an angular average and defining 

a(e) := (u(e)a(e)), 

Av{e) 



R(e) := u{e) - 

t 

„ ,Au(e). 



Using the distribution function a, introduced in the pre- 
vious section, we can write this expression as 



Sp n = -2eN(0) 

= -2e 2 iV(0)[$„-*„], 



e<&„ — / de ( u(fc, e)a(/c, e) 



(46) 



with the charge-imbalance potential determined by the 
formula 

poo 

# = (l/e) / de(u(k,e)a(k,e)). (47) 
Jo 

Relation of a(k, e) to the "clean limit" charge- 
imbalance distribution function ctk in k-space is given 
in Appendix D. The relaxation rate v(e) describes the 
relaxation into thermal equilibrium of the difference be- 
tween electron and hole- like quasi-particles due to impu- 
rity scattering. It replaces the relaxation rate l/r q in 
Ea. (jT2"|) . The r.h.s. of Eq.(g2J) contains the tunnelling 
of quasi-particle charge from the normal electrode into 
the superconducting layer due to the applied voltage. It 
is similar but not equal to the tunnelling current at the 
same energy e. In the tunnelling current (see below) the 
function R(e) is replaced by the density-of-states func- 
tion u(e). In the absence of impurity scattering one finds 
R(e) = l/u{e) in agreement with Refill 



(42) 



we obtain 

= WW)_ m£ _ v) _ p/(e + v) y (43) 
v{e) 

Finally from Eq.JSSJ) and Ea. (|4T)|) we obtain for the 
distribution function on the superconducting layer 

P(e) = \{P'{e-v)+l3'{e + v)l (44) 

where /3'(e) is the distribution function on the normal 
layer, which is assumed to be in equilibrium. This sim- 
ple relation is only true as long as we neglect inelastic 
scattering, which is necessary for a relaxation of quasi- 
particles into the condensate. 



E. Charge-imbalance 

Equation (|43|l is the main result of the paper. The 
function <5(e) describes the charge of quasi-particles with 
energy e on the superconducting layer. It is related to the 
quasi-particle potential \& introduced in the phenomeno- 
logical theory in the following way: The quasiclassical 
expression for the charge density is 



Sp-n 



-2eN{Q) 



e<5> r . 



V. CALCULATION OF THE TUNNELLING 
CURRENT 

In order to measure the charge-imbalance induced by 
the quasi-particle injection we also need an expression 
for the current, in particular between the normal contact 
and the first superconducting layer. Quite generally the 
current between neighboring layers n and n' in lowest 
order in the tunnelling matrix element can be written as 

Jn,n'(*)=+¥y] / ^ll^'fc'l 2 

(6%(k, t, ti)e" JT37 "."' (tl)/2 G< (fc, ti,t)e +iT3 ^- n + 1 ^/ 2 
G<(k, t, h)e- ir ^^y 2 G*(k, t u t)e+ iT ^'W 2 ) ^ 



+ c.c. 



(48) 



(45) 



In the stationary case (constant applied voltage) we ne- 
glect the time dependence of the Green functions on the 
central time and use a Fourier transformation with re- 
spect to the time difference as defined above. Further- 
more we restrict ourselves to the tunnelling between the 
normal layer and the first superconducting layer. Then 
only the normal Green functions contribute and we ob- 
tain (with v — £l n ,ri /2) the usual expression for the 
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quasi-particle current density 



'(«) 



-iV(O) 



(«i/ t (fc,fc / )ff R (*.e)s < (* / ,c-«)) 
+ („ t (fc, k')g<(k, e " «)») + c.c. (49) 



Here the prime denotes momenta and distribution func- 
tions on layer n' , and the double brackets denote an av- 
erage over both k and k' . 

Now we express the lesser functions g < by the Keldysh 
functions 



9< = l(9 K -9 R -9 A ) 



(50) 



and express the latter by the non-equilibrium distribution 
functions (3 and a, then 

2<= r +oc „ . „ 
3n,n<{v) = ~J-N{0) J de((u(k,k')u(k,e)u'(k',e-v) 

[[3'(e-v)-(3(e) + a'(e-v)-a(e)])) 

=-fnA V ) +in,n'(«). (51) 

This is another important result. The current between 
neighboring layers n and n' is the sum of two parts con- 
taining the distribution functions /3(e) and a(e) respec- 
tively. The current j°(v) is the quasi-particle current 
driven and created by the electro-chemical potential dif- 
ference between the two layers. The current j a (v) de- 
scribes the diffusion current driven by the charge imbal- 
ance. Both current contributions depend on the density 
of state u(e) of the two layers. 

For the further application to the tunnelling between 
a normal electrode and a superconducting layer it is con- 
venient to denote by j the current flow from the normal 
electrode to the superconducting layer and by V the volt- 
age drop in this direction (i.e. v = —eV). With u'(e) = 1, 
a' = 0, /3'(e) = /3 (e) = tanh(e/(fc B T)) for the normal 
layer, and exploiting the fact that u(e) and a(e) are even 
functions, while /3(e) is an odd function of e we obtain 
for the current contributions: 



f(V) = ^N{0)v t j de(u(k, e)) [/3 (e + eV) 



2e 



fa{e-eV)] =:jo(V), (52) 



J a (V) = --N(0)is t I de(u(k, e)a(k, e)) 
= : -o- V(V). 



(53) 



Here we have replaced the tunnelling rate by some av- 
erage v t = (v(k,k )) and have defined the ohmic resis- 
tance do = Ae 2 N (Q)v t / h. Then the current driven by 
the nonequilibrium distribution of quasi-particles can be 
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A=25 meV, v p =l meV 
T=4.2 K 

weak scatt. 

• strong scatt. 



10 20 
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40 



FIG. 2: Spectral function (tunnelling density of states) calcu- 
lated for a d-wave superconductor with impurity scattering. 



expressed by the quasi-particle potential 
1 



de(u(k, e)a(k, e))) 



(54) 



1 f v 

- / de^-(R(e))[0 o (e + eV) - /3 (e - eV)}. 
e Jo K e ) 



VI. RESULTS 

For our model calculations we use a d-wave order pa- 
rameter with an angular dependence of the usual form 
A(jfc) = A(T)(f)(k) with <p(k) = cos26» fe . The tun- 
nelling matrix element will be parametrised as v t {k,k') = 
Vi + f24>{k)<j){k'). This simplifies the averaging procedure: 
the parameter v\ enters the normal tunnelling probabil- 
ity while the parameter V2 determines the Josephson cou- 
pling. 



A. Tunneling current and charge imbalance 

In Fig. 2 we show well-known typical results for the 
spectral function (tunnelling density of states) u(e) at 
low temperatures with a self-consistently determined self- 
energy 17(e) for the two limiting cases of Born scatter- 
ing and in the unitary limit. Note that in the unitary 
limit the spectral function stays finite for e — > 0. This 
function will be needed as input for the following calcula- 
tions of the tunnelling current and the charge-imbalance 
relaxation. Fig. 3 shows the corresponding normalized 
tunnelling current jo (V) / oo calculated from Eq. (|52l be- 
tween a normal electrode and a superconducting layer 
as function of the voltage-drop V between the normal 
electrode and the superconducting layer. 
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FIG. 3: Tunnelling current jo/o"o[mV] as function of the ap- 
plied voltage. 
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FIG. 5: Charge imbalance potential generated on a supercon- 
ducting layer in contact with a normal electrode as function 
of the applied voltage (a) and the tunnelling current jo/o"o 
(b) 



FIG. 4: Relaxation rate for the charge-imbalance due to im- 
purity scattering as function of energy. 



voltage curve. 



Now we turn to the calculation of charge imbalance. In 
Fig. 4 the frequency dependence of the relaxation func- 
tion v(e) is shown for low temperatures. The frequency 
dependence reflects the available phase-space for elas- 
tic scattering processes from electron-like into hole-like 
quasi-particles at the energy e. With help of this function 
we calculate from Eq. (|54|l the charge imbalance potential 
'F generated on a superconducting layer in contact with 
a normal layer. Fig. 5a shows ^(V) as function of the 
voltage V between the normal electrode and the super- 
conducting layer. 

As V cannot easily be measured directly, we express 
it by the corresponding current density. The resulting 
function ip(j) defined by — *ff(V) with j — jo(V) 
given by Ea. (|52l) is shown in Fig. 5b. It depends on the 
ratio v t jv p between the average tunnelling rate and the 
potential scattering rate. The nonlinear dependence on 
the current reflects to some extent the nonlinear current 



B. Experiment with two normal electrodes 

Now we apply the theory to the basic experiment, 
where one superconducting layer is in contact with two 
normal electrodes. Through the first electrode with area 
F\ a current I\ is applied, which creates a charge imbal- 
ance ^> s on the superconducting layer. At a second elec- 
trode the voltage V% is measured with no current flowing. 
We assume that charge imbalance spreads evenly over the 
whole superconducting layer of size F . Then the charge 
imbalance potential created on the superconducting layer 
is given by $ s = {F\/ 'F)^{ji), where ji — I\jF\ is the 
current density through the electrode (1), with de- 
fined above and shown in Fig. 5b. 

In order to determine the voltage V-x measured at the 
second electrode we use the current equations (15^1 Q3J) 
and exploit the condition that no current is flowing, 
h = Jo(^2) — fo^s = 0, i.e. we have a compensation 
of the quasi-particle current driven by the voltage and 
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A=2S meV, v, / v p =0.01 
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FIG. 6: Voltage measured at the second electrode as function 
of the injection current through the first electrode. 



the quasi-particle diffusion current driven by the charge 
imbalance potential. Using for ^ s the value determined 
above we find 

de(u(k,e)[(3 (e+eV 2 )-(3 (e-eV 2 )} = -^V(ji)- (55) 

Thus we obtain V% as function of j\ . Results are shown in 
Fig. 6 for the case of a very small test electrode, F 2 <c F, 
Fx ~ F. Note that in this channel the electrochemical 
potential drop is only between the normal electrode and 
the first superconducting layer, all the other barriers in 
the stack have fl n , n +i = 0. Furthermore, if the elec- 
trodes are in equilibrium, the total voltage equals the 
total electrochemical potential, thus V% is the total mea- 
sured voltage. 



C. Double- mesa experiment 

In this experiment which is described in detail in Refi 
two small mesas are structured close to each other on top 
of a common base mesa and contacted with separate gold 
electrodes. Through the first mesa a variable current I\ 
is injected while at the second mesa the voltage V 2 is 
measured for fixed current I 2 . Normally the voltage V 2 
is independent of the current I\ as long as all junctions in 
the base mesa are in the superconducting state. In some 
cases, however, a small additional voltage AV 2 (I\) is ob- 
served. This happens, if the lowest junctions in the two 
mesas are in the resistive state (see Fig. 7b) and generate 
a charge-imbalance on the first common superconduct- 
ing layer of the base mesa, which then depends on both 
currents. 

If we want to apply the microscopic theory to this sit- 
uation we have to calculate the charge-imbalance on a 
superconducting layer in contact with another supercon- 
ducting layer with the barrier being in the resistive state. 
Furthermore we have to add to the tunnelling current the 



B 




FIG. 7: Position of resistive barriers in the double-mesa ex- 
periment 



average Josephson current. In systems with a large Mc- 
Cumber parameter this contribution is small. We will 
neglect it in the following. 

In a first approximation we may treat the junctions 
between the two mesas and the base mesa as if they were 
normal electrodes. In this case the only modification is 
the finite current through the second mesa. The charge 
imbalance potential generated on the first superconduct- 
ing layer of the base mesa is 



*b = 



T^(tt) + -ttVKtt) 



F Fx 



F^F 2 



(56) 



The current density through the second mesa, which is 
kept constant, is 



h = jo(V 2 ) - (Jq^b, 



(57) 



where jo(V) is the current-voltage function (|52|l . The 
voltage shift AV 2 (Ii) = V 2 (Ix) — V 2 , is the difference 
between the voltage measured at the second mesa for 
fixed current I 2 , when the last barriers in the first and 
second mesa are in the resistive state (Fig. 7b), and the 
constant voltage V 2 , when only the last barrier of the 
second mesa is in the resistive state (Fig. 7a). Expanding 
jo(V) around the voltage V 2 as jo{V 2 + AV 2 ) = jo{V 2 ) + 
AV 2 a(V 2 ) we obtain: 



(58) 



where the function ip(j) is shown in Fig. 5b. 

This approximation can be improved if we take into 
account that the last layers in the small mesas are in 
the superconducting state. Then we have to include the 
frequency dependent density of states u'(e) of these lay- 
ers both in the calculation of the quasi-particle current 
densities and in the charge imbalance function a(e). Let 
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FIG. 8: Modified charge imbalance potential generated on a 
superconducting layer by resistive junction as function of the 
tunnelling current. 



The corresponding voltage shift for a given current j 
is then 



A7 = 2 ^")» 



(61) 



where V is the voltage for a single junction. This shift 
can be measured directly as difference between the volt- 
age of two isolated junctions in the resistive state and 
two neighboring resistive junctions in the stack. This 
generalises our results obtained with the phenomenologi- 
cal theory! using an ohmic quasi-particle IV-curve. Note 
that this (very small) voltage shift in the IV-curves does 
not depend on the parameter a (called a in RefsF?1l2ll de- 
scribing charge fluctuations of the superconducting con- 
densate. These will be of importance in dynamical ef- 
fects like the Josephson plasma resonance and in optical 
experiments^ 



VII. SUMMARY 



us denote the modified quasi-particle current by jo(V) 
and the modified charge-imbalance potential function by 
ip(J) (explicite formulas are given in the Appendix), then 
the voltage shift AV^/i) measured at the second mesa 
is given by 



AV 2 (h 



(TO 



a(V 2 ) F V[ F 1 > 



(59) 



where cr(V) = djo(V)/dV). Its current dependence is 
proportional to the function -0(j) , which is shown in 
Fig. 8. The shape of the functions Fig. 5b and 
Fig. 8 is similar. The curve obtained in the weak scatter- 
ing limit (solid line) also has great similarity with the ex- 
perimental results. 6 The parameters used in our calcula- 
tions are typical for these materials. We did not perform 
a fit to the experimental data, since this would require 
an additional parameter for the Josephson coupling. Of 
course all the curves end at the critical current. 



D. Influence of charge-imbalance on the 
current-voltage curves 

The same formalism can be applied to calculate the in- 
fluence of charge-imbalance on the current-voltage curves 
in a stack of Josephson junctions. For instance, for a 
stack of junctions with one junction in the resistive state 
inside the stack (not adjoining the normal electrodes) the 
current is given by 



(60) 



where V is the voltage drop across this barrier. The fac- 
tor 2 comes from the charge-imbalance potential created 
on the two superconducting layers adjoining the barrier 
in the resistive state. 



In this paper we have developed a comprehensive mi- 
croscopic theory for stationary nonequilibrium effects 
in intrinsic Josephson systems starting from a non- 
equilibrium Green function theory for layered d-wave 
superconductors. We investigated the charge-imbalance 
generated on a superconducting layer by current injec- 
tion and derived results for the charge-imbalance dis- 
tribution function and the nonequilibrium quasi-particle 
current between superconducting layers. The theory 
uses basic nonequlibrium concepts developed earlier for 
bulk superconductors and is applied here to layered d- 
wave superconductors forming a stack of Josephson junc- 
tions. Specific for layered superconductors with small 
tunneling rate between the layers is the confinement of 
charge-imbalance on single superconducting layers. Spe- 
cific for d-wave superconductors with vanishing gap is 
the relaxtion of charge-imbalance due to elastic impu- 
rity scattering, which is the dominant relaxation mech- 
anism at low energies. In distinction to an earlier semi- 
phenomenological theory by the authors we considered 
here the energy dependence of the charge-imbalance dis- 
tribution function and its relaxation, leading to non- 
linear current-voltage realtions. 

We applied the theory to the calculation of nonequilib- 
rium effects in current injection experiments with 4 con- 
tacts. In particular, we calculated the voltage between 
a normal electrode and a superconducting layer for zero 
current as function of the current through a second elec- 
trode. This voltage measures directly the charge imbal- 
ance potential generated on the superconducting layer. 
We then applied the theory to recent double-mesa ex- 
periments. Thus we were able to explain the non-linear 
dependence of the voltage measured at one mesa on the 
current through the second mesa. 

The same formalism can also be applied to calculate 
the influence of charge imbalance on the current-voltage 
curves in a stack of Josephson junctions, which should 
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be observable as difference in the voltage between dif- 
ferent configurations of a given number of resistive junc- 
tions in the stack. We note that the shift of the chemical 
potential of the condensate leads do a redistribution of 
the voltage between different superconducting layers, but 
this has no influence on the total current-voltage curves. 
Charge oscillations in the condensate and the resulting 
coupling between the layers will however be important 
for dynamic effects like the dispersion of the longitudinal 
Josephson plasma resonance and in some optical exper- 
iments. These dynamic effects will be investigated in a 
forthcoming publication. 



functions we introduce a Fourier transform with respect 
to the time-difference r := t\ — t 2 keeping the central 
time t := (ti + t 2 )/2 fixed: 

G(k, t,e) = J drG(k, t+^,t- ^)e ieT . (A5) 

Finally we introduce the quasi-classical approximation 
by integrating over the energy = e k — fx keeping the 
direction k of the momentum fixed: 

g(k,t,e) = ( 9 f+ f ^j=^J d£ k G(k,t,e). (A6) 
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APPENDIX A: DEFINITION OF 
NON-EQUILIBRIUM GREEN FUNCTIONS FOR 
LAYERED SUPERCONDUCTORS 

Following the standard definitionsi£ii£ for non- 
equilibrium Green functions for superconductors we de- 
fine larger and lesser Green functions for layers n and nl 
(k = —k) in Nambu space by 

G> n ,(k,tv,k',h)= (Al) 

_ . / (Cnk<t(tl)cl, k ^(t 2 )) (Cnkt(t 1 )c n ,j., l (t 2 )) \ 

G< n ,(k,t 1] k',t 2 )= (A2) 
. ( (ci'fc'T^ 2 ) c " fc T(*i)) (Cn'fc'l^KfcT^l)) \ 

+l \-{cl kl ,{h)i H {h)) - w^fo)); 1 

from which we obtain the retarded advanced and Keldysh 
function by 

G R = Q{t 1 -t 2 ){& > -G<), 
G A = -e(t 2 -ti)(G> -G<), 
G K = G> + G<. 



(A3) 



For the different components in Nambu space the follow- 
ing notation is commonly used: 



G 



G F 
-F+ G 



(A4) 



For the average diagonal Green functions with n' = n 
only one index and one k-vector will be used. For these 



APPENDIX B: DERIVATION OF KINETIC 
EQUATIONS 

In our model the equations of motion for the Keldysh 
Green function on layer n in Nambu space are given by: 

i^T 3 G%(k, h,t 2 ) - H(k, h)G*{k, hM) 

- {t R G K + V K G A ] = 0, 

-i^-G^(k, t u t 2 )T 3 - G* (k, h,t 2 )H(k, t 2 ) 
ot 2 



{G £ + G £ } = 0, (Bl) 



where 



H{k,t) = -ecj> n {t)+i k -k k {t), 



(B2) 



with the clcctical scalar potential <fi n (t) , the charge of the 
electron (— e), and the order parameter matrix 



A fc (t) = 







A fc (t) 



-At(t) 



(B3) 



Here Ak{t) — Ake lXn ^ has a time-dependent phase. 
The constant amplitude of the order parameter is 
equal on each layer and has d-wave symmetry. In or- 
der to eliminate the time-dependent phase in the or- 
der parameter we make a gauge transformation c n k (t) = 
c n ,k(t) exp(ixn(t)). After this gauge transformation the 
new Green function G fulfills an equation of motion with 
H(k,t) = — e<& n (t) + 6c — Afc, the gauge invariant scalar 
potential <E> n (t) = 4> n {t) — hx n (t)/2e, and a real time- 
independent order parameter. 

The symbol {£G} denotes a convolution in time space 
of the self-energy and the Green function: 

{AB}(h,t 2 ) = J dtnA(t u H)B(t s ,t2). (B4) 

The self-energy, which will be discussed in detail later, 
contains random impurity scattering within the layer 
and tunnelling to the neighboring layers. The latter will 
treated in second order in the tunnelling matrix element. 
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A kinetic equation is obtained subtracting the two 
equations, introducing a central time t := (t\ -\-t%)/2 and 
taking the Fourier transform with respect to the time- 
difference r := t\ — t 2 : 

(B5) 



{± R G K + ± K G A - G R t K - G K ± A ). 



In this equation the kinetic energy drops out and we 
can perform the integration over keeping the direction 
k of the momentum fixed. We then obtain the kinetic 
equation for the quasi-classical Green functions 



. i d 

{ 2di +e)T39 



K 



i d 



e)g K r 3 



y 2dt 

- {h(k,t)g K - g K h(k,t)} 
{o R g K +o K g A -g R o K -g K o A } (B6) 



with 



h{k,i) = -e$„(t) - A(fc), 

A(fc)\ 
-A(fc) / 



A(fc) 



(B7) 



The curly brackets denote a convolution in time and fre- 
quency space: 



Note that all the contributions depending on the phase 
difference <p drop out. The remaining terms vanish if we 
neglect terms of order v t tt which are of higher order in 
the tunnelling probability. 

b) Tunnelling between two superconducting layers with 
the barrier in the resistive state. 

In this case also the anomalous Green functions con- 
tribute to the dc current. As this contribution involves 
the coherent part of the tunnelling matrix element and 
will be small in the limit of large McCumber parameters, 
we neglect it in the following. However, we have to take 
into account the frequency dependent density of states 
u(e),u'(e) on both superconducting layers. We want to 
apply the theory to calculate the charge-imbalance a on 
a superconducting layer n which is coupled on one side 
to a superconducting layer n' in a stack with barriers 
in the resistive state and on the other side to a super- 
conducting layer with the barrier in the superconducting 
state. As the layer n' is between two barriers in the resis- 
tive state with equal quasi-particle current we can neglect 
the charge-imbalance a' on this layer. Then we obtain 
the contribution to the scatterng term from tunnelling 
between n and n': 

III ~ 1% = -4*(^(«'(e -v)- u'(e + «))u(e)/3(e))' 
+4i{v t (u'(e - v)/3'(e - v) - u'{e + v)f3'(e + v))u(e))' 
-U{v t {u'(e -v) + u'(e + v))u(e)a{e))' , 



{AB}(t, e) = e*^ a « " 9 « a ' ^ 2 A(t, e)B{t, e). (B8) J*» - J% = M{v t v{e)) (u\e - v)(3'{e - v) 

-u'(e + v)(3'(e + v)))'. 



The equation l|B6|l is the starting point of our calcula- 
tions. 



APPENDIX C: COUPLING BETWEEN TWO 
SUPERCONDUCTING LAYERS 

a) Tunnelling between two superconducting layers with 
the barrier in the superconducting state 

Here we have £l n ,n' = but a finite constant phase 
difference <p := 7„ n /, and we obtain for the tunnelling 
self-energy 

of 0c, e)n, n i = - i(iy t g R {e)YT3 + sintp{v t f R (e))'n 

+ cos^ t /«( e ))'r 2 , (CI) 

of (k,e) = - 2^< ^ V(e) / 3 / ( e ))'T3 - 2i(v t u> '(e)a' '(e))' V 



(C4) 



For the charge- imbalance function a(e) on the supercon- 
ducting layer n we then find: 

2u p [u(e){u(e)a(e)) - {u(e))u(e)a(e)] 
Av(e) 



2v n 



u(e)a(e)) 



= -( tt (c) - jv(e))(u t (k, k') [u'(e - v)((3'(e -v)- (3(e)) 
+u '(e + v)(f3(e)-/3'(e + «))]>'. (C5) 



The formulas for the current densities obtained from 
Eq. (|51() now read in (symmetrised form) 



2h 



de(u(e)){u'(e + eV) 



(C6) 



2sm<p(v t v'(e)p'(e)yn 
2coscp(v t v'(e)(3'(e)yT 2 
2 cos tp(v t w' (e)ct (e)) 1 T\ 
2sm(p(v t w'(e)a'(e))'T 2 . 



.f(V) 



+ u'(e - eV))(f3 (e + eV) - (3 (e - eV)) =: j (V), 

f-\-QO 



-N(0)v t 



dea(e)(u'(e + eV) 



+ u'(e-eV)) =:-o *{V), 



(C7) 



(C2) 



For the contribution from tunnelling to the scattering 
term in the kinetic equation we then find: 



Tts jts 

jts jts 



u - ^2 - -8i(u t u'(e)u(e)a(e)y + Ai(v t u' (e)a' {t)u{e))\ 



with the charge imbalance function a(e) = (u(e)a(e)) 
given by 

fi(e) = Vt f^ {u'(e + eV) + u'(e - eV) 



2v(e) 



J% - 4{ = 8i(v t u>(e)a'(e)v(e))'. 



(C3) 



)(/3 (e + eV)-P (e-eV)). (C8) 
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The current equation for j a (V) defines a modified charge- 
imbalance potential 4 r (V), from which we obtain the 
function ip{j) = &{V) using the current-voltage relation 
j = jo(V). This function is shown in Fig. 8. 

APPENDIX D: CHARGE IMBALANCE OF 
QUASI-PARTICLES 

In order to make contact with the traditional theory 
of charge imbalance let us summarise the basic defini- 
tions and concepts of charge-imbalance as introduced by 
Tinkham and Clarka 10 ! 11 ! 31 ! 32 . 

In the BCS theory the total charge (in units of the 
electron charge, factor 2 from spin) is 

Q = 2J24 + («* - «*)/* =■ Q c + Q\ (Dl) 
k 

where u\ = |(1 + v 2 k = (1 + £k/E k )/2 are the 

usual coherence factors, = — M> ^fc = Vsl+^fc 
is the quasi-particle excitation energy. The first term is 
the condensate charge Q c , the second term is the quasi- 
particle charge Q* , fk is the quasi-particle distribution 
function. In equilibrium f k — l/{exp{E k /T) + 1) and 
hence Q* vanishes (for particle- hole symmetry). 

A noncquilibrium state can be described by a shift Sfi 
of the chemical potential and a change of the distribution 
function f k . A shift of the chemical potential leads to a 
shift of the excitation energy in k-space, — > e k — fi + 8fi 
and hence to a change of the condensate charge by: 

8Q C = 2N {0)8(1, (D2) 



where N{0) is the density of states for one spin direction 
at the Fermi energy. 

The quasi-particle charge can also be written as 

Q* = 2 Vkfk = - h< )' ( D3 ) 

k k 

where qu = £,k/E k can be considered as charge of a quasi- 
particle with momentum k. For each quasi-particle state 
with momentum k = k > and £ k> > (electron-like quasi- 
particle) exists a quasi-particle state with momentum 
k = fc< and £ k< = — £ k> < (hole-like quasi-particle) 
with the same excitation energy E k > and direction 
k. Thus Q* depends on the difference in the number of 
electron- and hole-like quasi-particles. While f k> — f k< 
describes charge-imbalance, the combination f k> + fk < , 
for instance, enters the self-consistency equation for the 
gap A fe . 

These distribution functions are introduced here for 
well-defined quasi-particles with infinite life-time. In a 
microscopic theory based on noncquilibrium Green func- 
tions these are replaced here by the frequency-dependent 
distribution functions a{e) and /3(e). In our case we find 
the following correspondence for e > 0: 

0{e = E k ) = l-f k> -f k< , (D4) 

a{e = E k ) = -q k {f k> - f k< ). (D5) 

Extended to the whole frequency range p(e) becomes an 
odd function, a{e) is an even function. 



* On leave from the Institute for Physics of Microstructures, 
RAS, Nizhny Novgorod, Russia 
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